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Abstract 

This paper is devoted to study some properties of the fc-dimensional Lyness' map 
F{xi, . . . ,Xk) = {x2, . ■ . ,Xk,{a + X]i=2 -^i)/^!)- main result presentes a rational 
vector field that gives a Lie symmetry for F. This vector field is used, for fc < 5, to 
give information about the nature of the invariant sets under F. When k is odd, we also 
present a new (as far as we know) first integral for FoF which allows to deduce in a very 
simple way several properties of the dynamical system generated by F. In particular for 
this case we prove that, except on a given codimension one algebraic set, none of the 
positive initial conditions can be a periodic point of odd period. 
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1 Introduction and main results 

Tlie second and more specially the third order Lyness' difference equations 

a + Vn+l , a + Vn+l + yn+2 . , , . „ 

yn+2 = and y„+3 = , with a > 0, 

Vn Vn 

have been considered as emblematic examples of integrable discrete systems, see for instance 
[T^ EH [23] • The dynamics of the above equations, or their associated maps, has been 



1 



the objective of recent intensive investigation. Nowadays the behaviour of the orbits is well 
known when positive initial conditions are considered, see [21 |5l El HOl El]- However few 
results have been obtained for negative initial conditions, see [6l [10] . 

On the other hand, the study of the 2 and 3-dimensional Lyness' maps is the starting 
point for the study of other birational integrable maps, see for instance O El [101 [131 [ISl [13 
[18] [23] . and [9] for a general paper on this topic. 

For k > 4, very few results, apart of the ones obtained recently by Bastien and Rogalsky 
jl], are known for the k-ih order Lyness' equation 

_ « + Li=l Vn+i ... 

Vn+k - • (ij 

Vn 

The main difference between the A; = 2 , 3 and the A: > 4 scenarios is that the first cases 
are integrable in the sense that the associated maps have 1 and 2 functionally independent 
first integrals, respectively (in this paper we say that a map F is integrable if it has k — 1 
functionally independent first integrals). It seems that this property is not shared for the 
Lyness' equations when fc > 4. 

Indeed, consider the fc-dimensional Lyness' map associated to the difference equation ([T|, 

F{xi, . . . , Xfc) = (x2, . . . , Xfc, " ^^=^ ^' I , with a>0, (2) 

which is a diffeomorphism from := {x = (xi, . . . , Xk) € M'^ : xi > 0, X2 > 0, . . . , > 0} 
into itself. It is well known that it has the following couple of functionally independent first 
integrals 

^i(x) =(^a + Y,Xi]{ Ylixi + 1) ] /{xi ■■■Xk) (3) 

and 

V'2(x) = ^a + ^Xj + xiXfcj + Xi + Xi+i) j /(xi • • -Xfc). (4) 

A third functionally independent first integral for k > 5 has recently been given in |11] . 
Moreover in that paper it is conjectured that for any k, the Lyness' map has up to E (^^) 
functionally independent first integrals, where £'(•) denotes the integer part function. The 
conjecture seems to be true for A; < 6, see again [IT] . 

The integrable structure for k = 2,3 implies that the dynamics of the maps studied in 
the above references is in fact one-dimensional. In any case, although for A; > 4 the above 
assertion seems not to be true, the existence of several first integrals reduces the dimension 
of the space where the dynamics takes place. If the above conjecture would be true, then 



i=l 
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k- E{^), would be, generically, the dimension of the invariant manifold given by the 
level sets of the first integrals. 

One geometrical object that has played a key role to understand the dynamics of a large 
class of 2 and 3-dimensional maps is the Lie symmetry of the map, [3 |8] . Recall, that a 
vector field X is said to be a Lie symmetry of a map G if it satisfies the condition 



X(G(x)) = (Z)G(x))X(x). 



(5) 



The vector field X is related with the dynamics of G in the following sense: G maps any 
orbit of the differential system determined by the vector field, to another orbit of this 
system, see [8] . In the integrable case, where the dynamics are in fact one dimensional, the 
existence of a Lie symmetry fully characterizes the dynamics. In [HI Thm. 1] it is proved 
that if G : ^ ^ ^ is a diffeomorphism having a Lie symmetry X, and such G preserves 7, a 
solution of the differential equation x = X(x), then the dynamics of G\^ is either conjugated 
to a rotation, conjugated to a translation of the line, or constant, according whether 7 is 
homeomorphic to S^, M, or a point, respectively. Other properties of the Lie symmetries of 
discrete systems are studied in |14| . 

One of the main results of this paper, is the following theorem where the expression of 
a Lie symmetry for the /c-dimensional Lyness' map is given. 

d 



Theorem 1. For /c > 3, the vector field X^ 
k-dimensional Lyness' map (flD, where 



.^^^Xi— — , is a Lie symmetry for the 

OXi 



Xi(x) 



(Xi + l) Ui=2i'^+Xi + Xi+i) (« + Ei=l 



Xi - X2Xk) 



n 



k 

i=2 



(6) 



Xi 



Xm(x) 



for all 2 < m < k — 1, and 



n 



(xfc + 1) Y\i=i{l + Xi + Xi+i) (a + EL2^« ~ ^i^fc-i) 



(7) 



(8) 



Xi 



Once we have the candidate X^ to be a Lie symmetry of the Lyness' map F the proof of 
Theorem [1] only will consist in checking that ([5| holds. We give here some hints of how we 
have found the above X^ . Observe that if there exists a vector field X^ , satisfying equation 
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^ for the fc-dimensional Lyness' map then 
/ Xi{F) \ / 10 



X2{F) 
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o \ / x,\ 

X2 



Xl 



Hence it is necessary that 

Xi+i=Xi{F), fori = l, 
and the "compatibility condition": 



Xl 



(9) 



Xk{F) 



1^1 + 



Xl 



.i=2 



Thus, the construction of a Lie symmetry is straightforward once the right expression of 
Xl, as a seed of equations ([9|, is obtained. In O [T] and [8l the expressions of X2 and X3 
are given. The idea of these papers is that these vector fields have to be multiples of VVi 
and VVi X VV2, respectively. These constructions are used to force F and to share the 
same set of first integrals. When /c > 4 we can not use anymore this idea because there are 
no enough first integrals. Nevertheless we observe that the first components of X2 and X3, 
are 



(xi + l)(a + Xl 



and 



{Xi + 1)(1 + X2 + X^){a + Xi + X2- X2X'i) 



X2 X2X3 

respectively. Thus it seems natural to try with Xi as the expression given in (|6| and, 
indeed, it works!. From this starting point, the proof for a given small k is only a matter of 
computations, while the proof for a general k is long and tedious, but straightforward. It 
is done in Section [2} We suggest to skip this section in a first reading of the paper. 

Another result that helps for understanding the dynamics generated by -F, when k is 
odd, is given in next proposition. In this result, the key point is the existence of a new (as 
far as we know) first integral for F^ = Fo F for any odd > 3. As we will see, our proof of 
the existence of this function is inspired in the paper |22J , where this first integral is given 
for A; = 3. 

Theorem 2. Set A; = 2^ + 1, x = (xi, . . . , X2t+i) and consider F from into itself. Then 
(a) The function 



W{-k) 



is a first integral of F'^ . 



ni=o(^2j+i + 1) 

ne ■ 
j=l X2j 



(10) 
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(b) For any £ > 2, the function := W + W{F), which is, 

11^=1 

is a first integral of F which is functionally independent with the first integrals V\ and 
V2 given in and Q), respectively. 

(c) The function W ■ W{F) coincides with the first integral of F, Vi given in (^. In other 
words, Vi = W ■ W(F). 

(d) The algebraic set ^ := {x G : VF(x) — W{F{'x)) = 0} is invariant by F. 

(e) If the map F has some periodic point of odd period then it has to be contained in Q. 

(f) Setting := {x e Q+ : ±(VF(x) - W{F{^))) > 0}, the map F sends G+ into Q- 
and viceversa, and both sets are invariant by F^. Furthermore, the dynamics of F^ 
on each of these sets are conjugated, being the map F itself the conjugation. 

(g) The measure 

mi{B) := / , , dx, 



Bn(x)(l^(x)-M^(F(x))) 
where n(x) = nf=i^«i invariant measure for F^, i.e. mi{F'^{B)) = mi{B), 

where B is any measurable set in Q^. 

(h) The measure 

m2{B) := / dx, 



JB n(x) 

is an invariant measure for F'^, i.e. m2{F'^{B)) = m2{B), where B is any measurable 
set in Q"*". 

We remark that, when k is odd, the first integral V3 given above coincides with the one 
given recently in [TT] . Observe also that the invariant algebraic surface G was already found 
in [7], but only for k = 3. Also, as we will see in Subsection 3.1 the function W is useful to 
make an explicit simple order reduction when we study the dynamics of F for k odd. 

Although, by using both theorems we have not been able to present a complete study 
of the higher dimensional Lyness' map, in next results we give some information about the 
invariant sets in the phase space when k = 4,5. We prove: 

Proposition 3. The vector field X4 given by equations (^^@ for k = 4, is a Lie symmetry 
for the 4-dimensional Lyness' map. Moreover, 'X./^iVi) = 0, for i = 1,2, and then the sets 
Ih,k '■= {^1 = ^} n {V2 = k} n Q'^ are invariant by F and by the flow 0/X4. 

Furthermore, if we assume that both first integrals intersect transversally on Ch,k, 0, 
connected component of I^k^ then C^k is diffeomorphic to a torus. 
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Proposition 4. The vector field X5 given by equations @^(^ for k = 5 is a Lie symmetry 
for the ^-dimensional Lyness' map. Moreover, X5(Vi) = 0, for i = 1,2,3, where 



and the sets Ih,k,i '■= {^1 = h} Ci {V2 = k} n {V3 = £} H Q+, which generically have at least 
two connected components, are invariant by F and by the flow ofX.^. 

Furthermore, if we assume that the three first integrals intersect transversally on Ch^k/? 
a connected component of Ih,k,e, then Ch,k,e is diffeomorphic to a (two-dimensional) torus. 

It is important to notice that in tlie above results we do not assert that most of the 
connected components Ch,k and Ch,k,i are tori, since we did not succeed to prove that over 
them the intersection of the energy levels at the whole sets C^^k and Ch^k/ are transversal. 
Hence it remains open to decide whether they are two-dimensional differentiable manifolds 
or not, and in the case that they are not differentiable manifolds to decide which are their 
topology. In any case, our result reduces the problem to a computational question. 



Figure 1. Projections into of the flow associated to the Lie symmetry X4, and the orbit 
of the Lyness' map, for k = 4 and a = 4, both with initial condition (1,2,3,4). 

Our numerical simulations seem to indicate that for A: = 4, all the generic level curves 
Ih,k, are connected. On the other hand, for k = 5, generically the sets Ih,k,i seem to have 
exactly two connected components. In Figures 1 and 2 we give a projection in of these 
surfaces. Indeed, in Figure 1 we represent, both an orbit of F and an orbit of X4 starting 
with the same initial condition, and in Figure 2 an orbit of F for k = 5. Notice that in 
both cases the behavior of the orbits seems to indicate that two, (respectively three) is the 
maximum number of independent first integrals for F when k = 4 (resp. /c = 5), as it is 
suggested in [TT] . 




+X2 X4 (1 + X2) (1 + Xi) {a + Xi + X2 + Xs + Xi + X5) 
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Figure 2. Projection into of the first 5000 iterates of the Lyness' map for k = 5 and 
a = 1, staring at (1,2,3,4,5). Odd and even iterates are in different connected components. 

A remarkable fact that Figure 1 shows is that, although for k = 4 the manifold Ih^k is 
invariant for both the map F and the flow of X4, the map F seems to send an orbit of X4 
to a different orbit of the vector field. Further numerical experiments seems to confirm this 
fact. Under this situation, we can not use the techniques developed in |8]. This fact makes 
more difficult the knowledge of the behavior of F restricted to each J/i ^ and is one of the 
important differences between the cases A; = 2,3 and k = 4. 



2 Proof of Theorem [1] 



As we have already explained in the introduction, the existence of a vector field = 
Yli=i-^i'£ir satisfying equation (5) for the /c-dimensional Lyness' map (2), is equivalent to 
the set of equations 

Xi+i = Xi{F), foii = l,...,k-l, (12) 

together with the compatibility condition: 

/ , v^fc \ , r fc 



Xk{F) 



_2 J^i 



Xl 



.i=2 



(13) 



The proof will consist in checking that the choice of X^ given in the statement satisfies 



equations ( 12 ) and ( 13 ). The result is straightforward for A; = 3, 4, 5 and we omit the details. 



So, from now on, we assume that k > 6. 
We proceed in two steps: 

First step: We will show that from the expression ^ of Xi as a seed of equations ( |12[ ) 
we obtain the expressions of X^i for m = 2, . . . ,A; — 1 and Xk given by equations ([7| and 
(|8]), respectively. 
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Second step: We will prove that the compatibility condition (13) is satisfied. 



First step: We start with some preliminary observations: 
Observation 1. Set Ki := Xi + 1 for i = 1, . . . , A; — 1. Then Ki{F) = Xj+i + 1. 
Observation 2. If we set Li := 1 + Xj + Xj+i, then for all 1 < i < k — 1, Li{F) 
1 + Xi+i + Xi+2 = Li+i, and 

r 1 , , 1 , , « + Ej=2^* a + Y!l^^Xi+xiXk 
Lk^i{F) = l + xk + Xk+i = l + Xk-\ — — = — . 

Xl Xl 

For this reason 
(a) If Ml := nti^(l + + Xi+i), then 

+ Yli=l Xi + XlXk) (ni=2^(l + Xi+l\ 

Mi{F) = ^ 



Xl 

(b) Setting := H; 

=i,iy^m-i,m(^ + + Xj-f i) for 2 < m < — 2, we obtain 

(fc-l \ k 

n (1 
i=2,i^m,m+l / i=l 

(c) If Mfc_i := nti^(l + + then Mfc„i(F) = 11^2 (1 + Xi + x^+i). 

fc-i 

Observation 3. Set A?" = a + ^ Xj — X2Xfc, then 

i=l 



N{F) = a + '^Xi - xsXk+i = a + '^Xi - xa ""^^^' 

j=2 

_ « + E,i2 



i=2 i=2 ^1 



Xl 



(Xl - Xs) = Xfc+l(xi - Xs). 



k 

Observation 4. Set i? = a + ^ Xj + xix^, then 

i=l 



= a + Xj + Xfc+1 + X2Xfc+i = (1 + Xl + X2) = XA;+i(l + Xl + X2) 

^-^ Xl 



«=2 



Observation 5. (a) For all 2 < i < — 2 set 5j = Xi_i — Xj+i, then Si{F) = Xj — Xi+2. 
(b) Set Sk-i = Xk-2 - Xk, then 

^ . , _ a + E^=2 ^» _ « + Ef=2 - XlXk-l 

Jk-n-T ) — Xk-i - Xk+i — Xk-i — . 

Xl Xl 



8 



If we now consider the seed Xi given by equation using Observations 1, 2a, and 3 
we obtain that 



(X2 + 1) lli=3{'^ + Xi + Xi+i) {a + Yli=lXi+XiXk){xi-Xs) 



Xi 



Now applying systematically Observations 1, 2b, 4 and 5a, we obtain that for 2 < m < 
k — 1 the component = Xm~i{F) is given by equation ([T]). 
Observe that in particular 



Xk-i = Xk^2{F) 



{xk-i + l) l\i=i{l + Xi + Xi+i) {a + Yli=iXi + xiXk){xk-2 - Xk) 



hence the term Lk~i = 1 + Xk-i + Xk does not appear and so, in order to compute Xk we 
need to use Observations 1, 2c, 4 and 5b, obtaining the expression of Xk given by (jsl). 



Second step (compatibility condition ( |13[ )). A simple computations shows that 

Xk{F) = ^ 



xi 



where 



a + y^^Xj 

\ i=l / 



k-1 
.i=2 



+ ^ Xj + xi a + ^ Xi - X2Xk 



1=2 



Another computation gives that 



Xi 



.1=2 



i=3 



B 



where 
B 

and 



(Xi + 1) ( a + Y.U2 ) \ti=2 ( a + T,i=l - ^2Xk ) + 



" + Y.i=i Xi + xiXk ) C - Xfc(xfc + 1) ( a + Y.i=2 ~ xiXk-i ) Yli^i L 



k-l 



C — ^ ^ Xm{Xm + ^)Sm^n 



m=2 



Recall that Lj = 1 + Xj + Xj+i for all i = 1, . . . , — 1, 5m = Xm-i — Xm+i and M„ 



Y[i=i ijLm-i + Xi + Xj+i) for 2 < m < A; — 1. Therefore we want to prove that A = B. 
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Step 2a. First we show that B contains L2 and L3 as a factors. To see this, it suffices 
to check that L2 and L3 are factors of C. Observe that if L2 = 0, then we can write 
L3 = — X2, and X3 = — (1 + X2) hence xi — X3 = cci + a;2 + 1. Therefore 

^l{L2=0} ~ Sm=2 ^m{Xm + ^)SmMfn 

= Em=2 XmiXm + l)(Xm-l - X^+l) (U.iZli^m-l,mi^ + + ^i+l)) 
= X2(l + X2){xi - X3)(X4 - X2) (ni=4 (1 + + Xi+ifj 

+X2il + X2)ix2 - X4)il +X1 + X2) (nf=i(l + + Xj+l)) = 0. 

So it is a factor of C. 

If L3 = 0, then L2 = X2 — X4, and X3 — X5 = —(1 + X4 + X5). Hence 

~ Sm=3 Xm{Xm + l)(x„j_i — X„j+i) ^f^j!^]^ „j(l + + Xj+i)^ 

= X4(l + X4)(X2 - X4)(l + Xi + X2) (^0^=4^1 + Xi + Xi^i)^ 
-X4(l + X4)(l + X4 + X5)(l + X2 + X-i){l + Xi + X2) X 

X (nt5(i+^i+^i+i)) = 0- 

Hence L3 is a factor of C. 

The above results imply that 

C = L2LsQk{xi,X2, X4, . . . , Xk), 

(observe that X3 does not appear in the expression of Qk)- Hence L2 and L3 are factors in 
the expression of B, and then 

B = L2L3 [-(xi + 1) (a + Ef=2 x^) [Ui=l L^i\ (a + Ef=i - X2Xk) 
Furthermore 

Qk = -2:2(1 + X2) (jliZl Li^ + ix2 - X4)Li (jliZl + X4(l + X4)Li {\[\zl + 

+ Em=^5 Xm{Xm + l){Xm-l - Xm+i)Li (ni=4J^m-l,m -^j) • 

Step 2b. Now we state the following claim, that will we proved at the end of the proof. 

Claim: For A; > 6, Qfc(a:i, X2, X4, . . . ,Xfc) = (^11^=4 -^*) [^ia;2^fc-i - Xk-iXkLi\. 



By using the claim, B = 11^=2 ^ 



Li=2 



D, where 



D = -(xi + 1) + Yli=2 ^i) + Yli=i - X2Xfc^ Lfc_i 



+ + Y.i=l Xi + XiXkj [xiX2-La;_i - Xfc_iXfcLi] 

-Xfc(xfc + l)[a + Y!i=2 Xi - xiXk-1^ Li. 
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= -LiXk 
= -LiXk 
= -LiXk 



Step 2c. Observe that 

-Xk{xk + l)(a + Yli=2 - ^i^k-ij Li - Xk-iXkLi (^a + Yli=i + xiXkj 
{xk + l)(a + Yli=2 - xiXk-ij - Xk-i (^a + Yli=i + xiXk^ 
Xk-i (a + Yli=i Xi + xiXk -xi- xiXkj +{l + Xk) [a + Y!i=2 ^i) 
Xk-l [a + Y!i=2 ^i) + (1 + ^k) [a + Ei=2 ^i) 
= -LiXkLk-i (^a + Yli=2 ^i) ■ 

Thus D also contains the factor L^-i, and therefore B 

E = [-(Xi + l)(a + Yli=2 Xi) (« + E^t"l^ - X2Xk) 

+ (a + Yli=i Xi + xiXk) xiX2 - XkLi {a + Y!i=2 



E, where 



Step 2d. Is not difficult to check that E is a quadratic polynomial in X3. Now we see that 
E van 
when 

/ k-l 
.tq = .xq o := — I a A- > .r» -I 

1 + Xi 



E vanishes either when xs = xs^i := —a — X^^^^ Xi (that is when a + X^^L^ Xi = 0), and 

/ , , 3^2(1 - xiXk)\ 

X3 = X3^2 ■= - \a + 2_^Xi+ ' 



\ i=4 
In the first case we have 

^'{a+E^-i ^i=0} ^ ^1(^1 + ^)i-Xk - X2Xk) + x\x2Xk + LiXlXk 

= xiXk [—{xi + l)(a;2 + 1) + X1X2 + 1 + xi + X2] 

= xiXk [—1 — xi — X2 — X1X2 + X1X2 + 1 + xi + 2:2] = 0. 



In the second case, since a + 5^j=3 Xi = X2{x\Xk — l)/{x-i + 1), we have 

[xiXk - 1) V ^ ^\( X2{xiXk - 



T^l / . X2{xiXk - l)\ , (x2{xiXk-l) ^ 

^l{c.3=X3.2} = - ^2 + (^1 + 1) + Xl+X2-Xk- X2Xk 



( X2{x\Xk — 1) \ 

-\-XxX2 Xl + X2 H — V XlXk 

V Xl + 1 ) 

( X2(xiXfe-l)\ 1 

-LiXk X2 + — ; =: — j- F, 

\ + 1 / Xl + 1 

where (after some computations) 

F = -X2a;i(xfe + l)Li(xi - Xfe) + xfx2-Li(xfe + 1) - -LixiX2Xfe(xfe + 1) 

= XlX2Ll(Xfe + 1)[-(X1 - Xk) + Xl - Xk] = 0. 

Therefore E|{,3^,3 = 0. 

In summary, as a quadratic polynomial in X3, E factorizes as 

E = -(xi + 1) (x3 + a + Ejii,i^3 Xi) (x3 + a + J^fj^ Xi + (2:2(1 - a;iXfc))/(xi + 1)) 
= - (a + Xi) Ux3 + a + YliZl Xi) (xi + 1) + X2(l - a;iXfc)^ 
= - (a + Yli=i Xi) {{xi + 1) (a + Yli=3 Xi) +X2- xiX2Xk 



= - (a + E,^! Xi) (a + ^^^2 Xi + Xl (a + ^^^3 
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So, finally, we get that 



B 



as we wanted to show. 

To end the proof it only remains to prove the claim. We proceed by induction. That it 
is true when /c = 6 is straightforward. Assume now that the claim is true for Q^, then 

Qk+1 = -X2{1 + X2) \ U'i=iLi \ + (x2 -X4)Li [W^^^Li] +X4(1 + X4)Li (nLs^O + 



+ X]m=5 Xm{^ + Xm){Xm-l Xm+l)Li \ L, 

= LkQk + Li {jX^Zl Li) (xfc(l + Xfc)(xfc^i - Xk+l)) ■ 
By using the hypothesis of induction, we obtain that 
Qk+1 = Lk f Li) [xiX2Lk^i - Xk-iXkLi] + 

Ll (ni=4 Lij {Xk{l + Xk){xk-1 - Xk+l)) = 

(jli=i Li) {xiX2LkLk-i + Li [xk{l + Xk){xk-i - Xk+l) - Xk^iXkLk]) ■ 

Thus 

Qk+1 = Lk {j\^Zl Li) [xiX2Lk-i - Xk-iXkLi] + Li {\{\zl Li) {xk{l + Xk){xk-i - Xk+i)) 
= (l\i=i Li) {xiX2LkLk-i + Li [xk{l + Xk){xk-i - Xk+l) - Xk-iXkLk]) ■ 

An easy computation shows that Xk{l + Xk){xk~i — Xk+i) — Xk-iXkLk = —Lk-iXkXk+i, 
and the result follows. Therefore, Theorem [T] is proved. 

3 Geometrical issues in the odd case. 

Before proving Theorem |2] and Propositions [3] and |4] we need some preliminary results. 

v) 

Recall that a map H which is a first integral for := G o • • • o G is also called sometimes 



a p— first integral^ or simply, for short, a p— integral, of G, see [23j. 

The following lemma, which is very easy to prove, gives light on one utility of p— first 
integrals, specially if they are not symmetric functions of their arguments. 

Lemma 5. Let H be a p-integral of a map G. Then then for any symmetric function of p 
variables S, the function Vs := S{H, H o G, H o G"^, . . . ,H o G^^^) is a first integral of G. 

Lemma 6. Set k = 2i + 1, x = {xi,X2, ■ ■ ■ ,Xk) G and let W be the function given 



in (10). Then, if we define the polynomial 



^(x): ={Y{tl^^)[W{^)-W{F{x))] = 

= W]=0 X2j+l{x2j+l + l)-{a + Y.YJ{^ X^ n5=l X2j{x2j + 1), 

it holds that Z{F{x)) = det{DF{x))Z{x). 
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Proof. In Theorem |2] (a) we will prove that W is a 2-integral of F. Thus, if we define 
Z ■=W-W{F), then 

Z(F(x)) = -Z(x). (14) 

From the above equality it is clear that {x G : Z(x.) = 0} = {x G Q"*" : Z(x.) = 0} is 
an invariant hypersurface by F. Notice also that if we define n(x) := 11^=1 ^j' holds that 

n(F(x)) = " + ^^=2 ^^ n(x) = - det(DF(x)) n(x). (15) 



Since ^(x) = n(x)Z(x), by using equalities (14) and (15), we obtain that 



Z(F(x)) = n(F(x))Z(F(x)) = det(DF(x))n(x)Z(x) = det(L>F(x))Z(x), 
as we wanted to see. | 

Proof of Theorem [2| (a) The proof of the equality VF(F^(x)) = VF(x) is straightforward. 



We have obtained the expression (10) inspired in the results of [22 . In that paper it is 
proved that 

{Vn + l)(j/n+2 + 1) ^ {yn+2 + l)(yn+4 + 1) 
Un+l yn+3 

where {yn} is the sequence given by the 3-rd order Lyness' recurrence yn+3 = (a + yn+i + 
yn+2)/yn- Notice that this property is equivalent to say that for A; = 3, is a 2-integral 
for F. 

(b-c) By applying Lemma [s] with S{u,v) = u + v and S{u,v) = uv, we obtain the first 
integrals V3 and Vi, respectively. The functionally independence of Vi, V2 and V3, for ^ > 2, 
follows from straightforward computations and it is already established in jllj . 

(d-f) From Lemma [6] we know that Z{F{'x.)) = det(I?F(x))Z(x), where recall that 
^(x) = iULi^iWi^) - ^(^(x))]. Note also that 



det(DF(x)) = (-1) 



j.a + X2-\ h Xk-i 

4 



Since when k is odd det{DF) < on Q"*", equation Z{F) = det{DF)Z means that 

g = {x G Q+ : iy(x) = W{F{x))} = {x G Q+ : Z(x) = 0} 

is invariant by F and that F maps the region {x : Z{x) > 0} into the region {x : Z{x.) < 0} 
and viceversa. Furthermore it implies that the dynamics of F'^ on each of these sets are 
conjugated, being the map F itself the conjugation. Moreover, any periodic orbit with odd 
period must lie in ^, as we wanted to see. 
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(g-h) By using the Change of Variables Theorem it is easy to see that if G is a dif- 
feomorphism of li^ and on this region /x is a positive function that satisfies /i(G(x)) = 
det(L'G(x))/i(x), then 

[ 1 

m(B) = / ——— dx 
Jb 

is an invariant measure for G. By Lemma [6]we know that Z(F(x)) = det(DF(x))Z(x) and 
by equality (15), that n(F(x)) = — det(Z)F(x))n(x). By using these results we have that 
Z(F2(x)) = det(L>F2(x))Z(x) and U{F'^{x)) = det{D {x))U{x) , being both equalities in 
the corresponding domains, which are invariant by F'^. Hence (f) and (g) follow. | 



3.1 Order reduction 

Recall that in the previous section we have seen that when k = 2i + 1, the regions {x : 
Z(x) > 0} and {x : Z{'x.) < 0} are invariant by -F^, and the dynamics on both region 
are conjugated. This observation allows us to give a new application of the invariant W. 
We can reduce the study of the dynamics of F on {x : Z{x.) ^ 0} to the study of a new 
(fc — l)-dimensional map, having one more parameter. 

For instance for n = 3, we get that W{x,y,z) = (x + l)(z + l)/y, and hence any 
admissible level surface {x : W{x.) = w}, w ^ 0, can be described as y = k{x + l)(z + 1), 
where k = l/w. Therefore 

^21 / ^ a + z + k{x + l){z + l) a + k + z{k + l] 



l{l^=-}^'^'^'^^ " V ' X ' kx{z + l) 

and we can reduce the study of the dynamics of to the study of the reduced map 

~, , / a + k + z{k + iy 
kxiz + 1) 

or, equivalently, the study of the second order difference equation 

a + k + yn+i{k + 1) 

yn+2 = — 7—f —T\ — ' 

kyn[yn+i + 1) 

as in [22' Eqs. (8) and (9)]. The dynamics of this difference equation is studied in [8l Ex. 
3]. Equation (16) is sometimes called generalized Lyness' recurrence, see |20j . 

For n = 5, the integral of F"^ is W{x,y, z,t, s) = {x + l)(z + l)(s + l)/{yt). Again any 
admissible level surface {x : W{x.) = w}, w ^ 0, can be described as t = k{x + l){z + 
l)(s + l)/y, where k = l/w. Therefore proceeding as before we can reduce the study of the 
dynamics of to the study of the reduced map 

~ / p2{z,s;k)x'^ +pi{y,z,s;a,k)x + po{y,z,s;a,k) 

F2{x, y, z,s)= [ y, z, t, 2 

\ xy' 
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where P2{z, s; k) = k {s + 1) (z + 1), pi{y, z,s;a,k) = 2 k {s + 1) (z + 1) + y {a + s + z), and 
Po{y, z,s;a,k) = k {s + 1) {z + 1) + y {z + s + a + y), or equivalently to the difference equa- 
tion 

P2{yn+2,yn+3; k)yn^ + Pliyn+1, yn+2, Vn+S', a, k)yn + Po{yn+l, yn+2, yn+S] a, k) 
yn+i = 2 • 

Clearly, the described procedure can be generalized to higher dimensions. 



4 Dynamics of the low dimensional cases 

This subsection is devoted to prove Propositions |3] and |4| 

Along the section we will use a straightforward consequence of the following result: 

Theorem 7 (Bastien and Rogalsky, Let x be the fixed point of F in Q+. For any 
h > Vi(x), the level sets {Vi = h} n are homeomorphic to S'^^^. 

Corollary 8. Let K ^ ^ he by the intersection of some level sets of different first integrals 
of F, including V\ among them. Then K n is a compact set, invariant by F. 

For the sake of completeness and to compare with the cases /c = 4, 5, we start by recalling 
some results for the case k = 3, most of them already proved in [7]. 

4.1 The 3-dimensional map 

For A; = 3, 

Z7/ ^ f a + y + z 
F{x,y,z) = {y,z, 



X 

the Lie symmetry given in Theorem [T] is 



X 



3 



x{x + + y + z){a + X + y - yz) ^ + y{y + l){x - z){a + X + y + z + xz) ^ 
+z{z + 1)(1 + X + y){a + y + z-xy)-§^] /{xyz) 



and since lL^{Vi) = 0, for i = 1,2, the functions Vi and V2 given in ([s]) and Q are first 
integrals for X3 and F. Also, by Theorem [2] W{x, y, z) := (x + 1) (z + 1) /y is a 2-integral 
of -F; ^ = {x G : Z(x.) = 0} is invariant by F, where 

Z(x) = x{x + l)z{z -|-l)-(a-|-a;-|-y-|- z)y{y + 1); 

and F maps := {x G : Z(x.) > 0} into := {x G : Z{x.) < 0}, and viceversa. 

Let X be the fixed point in Q+, of F. Set h > Vi(x), k > V2(x). In [7^ it is proved that 
{Vi = h} n and {V2 = k} f] are diffeomorphic to two dimensional spheres, and their 
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transversal intersections are formed by exactly two disjoints curves, both diffeomorphic to 
circles. Their non transversal intersections correspond to: 

(a) The 2-periodic points of F (which are equilibrium points of X3) given by the curve 
C := {{x, {x + a)/{x — 1), x) \ X > 1}. 

(b) The levels placed at G. Those placed at ^ \ {x} are formed by exactly one curve, 
diffeomorphic to a circle. 

Finally, since X3 is also a Lie symmetry of F^, as a consequence of [71 Thm. 18] or [HI 
Thm. 1], we know that the map restricted to each of the sets {Vi = h}n{V2 = k}, which 
is not simply a fixed point of F^, is conjugated to a rotation. Further discussion about the 
possible rotation numbers can be found in [7]. 

4.2 The 4-dimensional map. 

Proof of Proposition |3j Equations ([6])-([8]) give the following Lie symmetry for the 
4-dimensional Lyness' map: 

X4 = [x {x + 1) {1 + y + z) {1 + z + t) {a + X + y + z - yt)£ 
+y {y + 1) {1 + z + t) {a + X + y + z + t + xt) {x - z)^ 
+z {z + 1) {1 + X + y) {a + X + y + z + t + xt) {y - t)£ 
-t {t + 1) {1 + X + y) {1 + y + z) {a + y + z + t - xz)§j] /{xyzt). 

A straightforward computation shows that the Lie symmetry satisfies X4(Vi) = 0, for 
i = 1,2. Hence, the orbits of both F and X4 generically lie in a two dimensional surface of 
the form Ih,k ■= {Vi = h} n {V2 = k} n Q+. 

Let Ch,k be a connected component of Ih,k- From Corollary |8] we know that Ch.k is 
compact. If {Vi = h} and {V2 = k} intersect transversally on C^^k: then for all points 



v^i A VV2 = muv2)y - {Vi)y{V2U£A ^ + [{v,uv2)z - {V^),{V2U£ A £ 

+[{VMV2)t - {VMV2).]i, A I + m)y{V2)z - {V,),{V2)y]l^ A I 
+ [{Vi)y{V2)t - {VMV2)y]i: A I + m),{V2)t - {Vi)t{V2)z]§-^ A |, 




{Vi)^ {Vi)y {Vi), {Vi)t 

{V2), {V2)y {V2)t 




This fact implies that the dual 2-form associated to the 2-field 



given by 



m),{V2)t - {VMV2)z]dxdy + [{Vi),{V2). - {ViUV2)z]dxdz 
+ [iVi)ziV2)y - iVi)y{V2)z]dxdt + [iVi)tiV2)^ - {Vi),{V2)t]dydz 
+ [iVMV2)y - iVi)y{V2)t]dydt + [iVi)M)y " {Vi)y{V2Udzdt 
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is nonzero at every point of C/^^fc, and therefore it is orientable, see [H Sec. 2.5] or [21] Sec. 
X, §1.]. 

Let X^''^ denote the restriction of the vector field X4 to the invariant surface C^^k- Some 
computations show that the unique equilibrium point of X4 in Q"*" is the fixed point of F. 
Hence X^''^ has no equilibrium points, and therefore the Poincare-Hopf formula gives: 

= i(X^''=) = x{Ch,k) =2-2g 

where i(X^'^) denotes the sum of the indices of the equilibrium points of X^''^ in Ch,k and 
xiCh,k) and g are the Euler characteristic and the genus of the surface Ch,k, respectively. 
Hence the genus of Ch,k is one. An orientable, compact, connected surface of genus one is 
a torus, as we wanted to prove. | 



4.3 The 5-dimensional map 

Proof of Proposition |4[ By Theorem |2] we know that 

r,r, , , (X+1)(Z + 1)(S+1) 

W{x,y,z,t,s) := , 

yt 

is a 2-integral of F. Moreover, the 3 functionally independent first integrals of F given in 
(§, (H) and ([ll|, are 

{a + x + y + z + t + s){x + l){y + l){z + l){t + l){s + l) 

Vi{x,y,z,t,s) = , 

xyzts 

, , {a + x + y + z + t + s + xs){l + x + y){l + y + z){l + z + t){l + t + s) 

V2{x,y,z,t,s) = , 

xyzts 

, , x{x + l)z{z + l)s{s + l) + {a + x + y + z + t + s)y{y + l)t{t + l) 

V3{x,y,z,t,s) = . 

xyzts 

Some tedious computations show that the hypersurface G is in the locus of non-trans- 
versality of the three level sets of the integrals Vi, i = 1,2,3 in Q^. Recall that precisely, 
G = {x £ Q'^ : Z(x.) = 0} is invariant by F, where 

Z(x) = x{x + l)z{z + l)s(s + l)-{a + x + y + z + t + s)y{y + l)t{t + 1), 

and that F maps g+ = {x G Q+ : Z(x) > 0} into = {x G Q+ : Z(x) < 0}, and 
viceversa. 

Equations ([6])-([8]) give the following Lie symmetry for the 5-dimensional Lyness' map: 
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a 



X5 = [x[x + l){l + y + z){l + z + t){l + t + s)[a + x + y + z + t-ys)Q^ 

dy 



+y{y + l){l + t + s){l + z + t){a + x + y + z + t + s + xs){x-z) 



+z {z + 1) {1 + X + y) {1 + t + s) {a + X + y + z + t + s + xs) {y - t) ^ 
+t{t + l){l + x + y){l + y + z){a + x + y + z + t + s + xs){z-s)f^ 
-s{s + l){l + x + y){l + y + z){l + z + t){a + y + z + t + s-tx)§^] /{xyzts). 

Again, direct computations show that Ji-^iVi) = 0, for i = 1,2,3. Hence the orbits of 
both F, and X5 he in a two dimensional surface of the form Ih,k,i '■= {^1 = h} Ci {V2 = 

k}n{V3 = i}nQ+. 

Let Ch,k,i be a connected component of Ih,k,i- From Corollary [8] we know that Ch^k/ is 
compact. If {Vi = h}, {V2 = k} and {V3 = k} intersect transversally on Ch,k,ii then for all 
points 

^ {Vl)x {Vl)y {Vl), {Vl)t {Vl)s ^ 

Rank {V2)x {V2)y (1^2). {V2)t {V2)s =3. 

V (^3)x {V3)y {V3)t (^3). y 

Similarly than in the case /c = 4, this fact implies that the dual 2-form associated to the 
3-field VVi A W2 A VV3 is nonzero at every point of Ch^k/-, and therefore this set is a 
two-dimensional orient able manifold. 

It is not difficult to check that all the equilibrium points of X5 in are the points of 
the curve 

f f 2x + a 2x + a \ 

£ = < X = X, , X, , X with X > 2 

[ \ X — 2 X — 2 / 

which contains a continuum of two periodic points and the fixed point of F. Moreover C 

belongs to the locus of non-transversality of the integrals Vi, and V2 in . 

The above observation implies that X5 has no equilibrium points in any level set Ih.k/ 

where the three first integrals intersect transversally. Therefore the Poincare-Hopf formula 

gives 

= i(xf = x(c,,fc,,) = 2-25. 

for each connected component Ch,k,e of such a level set (where Xg''^'^ is the restriction of 
X5 to Ch,k,e)- 

Hence g = 1, which implies that Ch.k/ is a torus (since it is two-dimensional, orientable, 
compact, manifold of genus one), as we wanted to proof. 

Finally, observe that the fact that F maps into Q~ , and viceversa, implies that most 
Ih,k,t have at least one connected component on each set. In fact, it seems that each Ih,k,ij 
not included in F, has exactly two connected components in Q"*", as it happens when A; = 3. 
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See Figures 2 and 3 for two illustrations of this assertion. 




Figure 3. Projection into of the first 10^ iterates of the Lyness' map, for k = 5 and 
o = 4, starting at (1, 2, 3, 4, 5). Odd and even iterates are in different connected components. 

From the above proof it is clear that if some Ih,k,e cuts C then the first three integrals 
do not cut transversally on it. Let us see that in general Ih,k,£ H £ = 0. This will be a 
consequence of the shape of the function 

, , / 2x + a 2x + a \ 

vi{x) := Vi X, —,x, —,x , X > 2, a > 0. 

\ X — 2 x — 2 J 

This function has a global minimum at the coordinate given by the fixed point x = 2 + 
\/4 + a, and lim vi{x) = lim vi{x) = +oo. 

Thus, given h > vi{2 + \/4 + a) there are only two solutions xi{h) < 2 + \/4 + a < X2{h) 
of the equation vi{x) = h. Now set 

, N f 2x + a 2x + a \ 

Vi{x) :=Vi [x, — — zro~'^l' x>2,a>0, 

for i = 2,3. Then, for any value of k and I satisfying k {i'2(a;i(/i)), f2(x2(/i))} oi i ^ 
{f3(a;i(/i)), U3(x2(/i))}, we obtain that Ih.k.t-, does not intersect C. 



5 Conclusions 

Several properties for the fc-dimensional Lyness' map F have been given, like the existence 
of a Lie symmetry for F and of a new and simple invariant for F^. This Lie symmetry 
together with the new invariant give information for k = A and 5 about the dynamics and 
the topology of the level surfaces where the dynamics of F is confined. Some general results 
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for k odd have been also presented. However, on the contrary that happens in the cases 
k = 2 and 3, the numerical explorations indicate that for k = 4,5 the orbits of the map are 
not contained in the orbits of the flow of the Lie symmetry with the same initial condition, 
although they are placed in the same manifold, which has dimension k — E{^^^). This is 
an obstruction to apply the theoretical tools developed in [S] . 

Numerical simulations seem to show that for some initial conditions the projection in 

of the iterates of F when A; = 6, 7 fill densely a 2-dimensional manifold, indicating that 
probably they live in a 3-dimensional manifold of and M'^. These facts are coherent with 
the conjecture of [llj about the number of independent first integrals of the Lyness' maps, 
and show that for k > 6 the dynamics are much more complicated. 

When k = 2£, the simplest scenario that we imagine for the dynamics of the /c-dimensional 
Lyness' map is that most of the orbits lie on invariant manifolds which are diffeomorphic to 
^-dimensional tori, S x ■ ■ ■ x S . On the other hand, when k = 2i + 1, most of them lie on 
two diffeomorphic copies S x ■ ■ ■ x S , separated by the invariant set Q. Moreover these 
orbits jump from one of these tori to the other one and viceversa. 

In any case, much more research must be done in order to have a total understanding 
of the dynamics and the geometrical structure of high dimensional Lyness' maps. 
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